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Abstract 

In this paper we present an integro-differential diffusion equation for continuous time random 

walk that is valid for a generic waiting time probability density function. Using this equation we 

also study diffusion behaviors for a couple of specific waiting time probability density functions such 
t— » - 

as exponential, and a combination of power law and generalized Mittag-Leffler function. We show 
f^j \ that for the case of the exponential waiting time probability density function a normal diffusion 

is generated and the probability density function is Gaussian distribution. In the case of the 
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combination of a power-law and generalized Mittag-Leffler waiting probability density function we 
obtain the subdiffusive behavior for all the time regions from small to large times, and probability 

qq \ density function is non-Gaussian distribution 
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I. INTRODUCTION 



Diffusion is a ubiquitous phenomenon, and it is one of the fundamental mechanisms 
for transport of materials in physical, chemical and biological systems. The well-known 
example of a diffusion process is the Brownian motion. Diffusion processes are classified 
according to their mean-square displacements. The normal diffusion is that the mean- 
squared displacement grows linearly with time, and in other situations the processes are 
said to exhibit anomalous diffusion. Nowadays, there are several approaches to describe 
anomalous diffusion processes, and they can be applied to many situations of natural systems 



9). One of the most interesting features incorporated into these approaches is the memory 



effect. In particular, the memory effect incorporated into the Langevin approach, referred to 
as generalized Lan gev in equation (GLE) |9|, can be associated to the retardation of friction 



and fractal media 



10 



111 ]. Moreover, according to the fluctuation-dissipation theorem 



Q. 



the internal friction is directly related to the correlation function of the random force. 

In many situations, a finite correlated noise is necessary for describing the real systems 
in equilibrium states. However, in order to describe anomalous diffusion, a nonlocal friction 
should be employed so that it satisfies the fluctuation-dissipation theorem. For instance, 
anomalous diffusion processes have been observed in a variety of systems such as bacteria 



cytoplasm motion 12], conformational fluctuations within a single protein molecule 13] 
and fluorescence intermittency in single enzymes [lj]. These processes have been described 
by the GLE, and the memory effect has also been shown in the generalized Fokker-Planck 
equations (GFPE) |9|. 

The continuous time random walk (CTRW) of Montroll and Weiss [15] has been employed 



to describe anomalous diffusion 
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221 ] . The CTRW with a power-law waiting time proba- 
bility density function (pdf ) 23|, |24| was also linked to the following fractional Fokker-Planck 
equation, 



where 



dp(x,t) 

dt 



Dl- a p(x,t) 



\l-a 



d 2 



-nDr a K a —p(x,t), 



1 d 



p(x,t X 



A-a 



dt x 



(1) 



(2) 



T(a)dtJ (t-h)' 
is the Riemann-Liouville fractional derivative and V(z) is the Gamma function. p(x,t)dx is 



the probability for finding particle in position between x and x + dx at time t. 



The CTRW model may be described by a set o 



appropriate generalized master equation 

n 

equation in Fourier-Laplace space [4|: 
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Langevin equations 
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261 ] or an 



281 ] . The pdf p(x,t) obeys the following 



p{k,s) 



l ~9(s) p {k) 



(3) 



s 1 — ip(k, s) 

where p Q (k) is the Fourier transform of the initial condition p (x), ij){x,t) is the jump pdf 
and g(t) is the waiting time pdf defined by 



g(t)= / ^(x,t)dx. ( 4 ) 

J— oo 

Moreover, from the jump pdf we also have the jump length pdf defined by 

POO 

cj)(x) = / ij;(x,t)dt . (5) 



The continuous time random walk model can be simplified through the decoupled jump 
pdf ip(k, s) = (p(k)g(s). Furthermore, jump length pdf can be approximated as following [4|: 



<j)(k) ~ 1 - Dk 2 + 0(k 4 



(6) 



Different types of the CTRW models are specified through specifying the wai ting time pdf 



The CTRW model is also connected to a class of Fokker-Planck equations [4, |29j. Further- 
more, in the CTRW model, solutions for p(x, t) under the condition of long-tailed waiting 
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32|. 



time pdf can be found in Refs. 

Substituting Eq.fjH]) into Eq.(j3]), we have 



p(k,s) 



Po(k) 



(7) 



s 1 - (1 - Dk 2 ) g(s) ' 

n 

When Cakir et al. [33] studied trajectory and density memory, the Eq.(43) in their paper 
is same as our Eq.flT]) with consideration of p {k) — 1. In their study, the initial condition 
is that when t = 0, particle is at origin, and therefore p (k) = 1. However, Eq. f^jty can be 



applied to cases without requirement of p (k) = 1 that appeared in Cakir et al.' paper 33 ] . 
In particular, the CTRW model can be classified by the characteristic waiting time T and 
the jump length variance S 2 defined by 



T 



tg(t)dt 



(8) 



and 

x 2 (p(x)dx . (9) 



oo 
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n 

For finite T and E , the long-time limit corresponds to the Brownian motion [4]. Although 
Eq. (j7j) is valid for a finite jump length variance, anomalous diffusion can be produced by 
Eq.(CJ) with appropriate choices of waiting time pdf. 

The aim of this work is: (1) from Eq.(J7]) to investigate general behavior of particle's 
diffusion without any specified waiting time pdf; (2) study the effects of different waiting 
time pdf 's on the behavior of particle's diffusion in the framework of CTRW. This paper is 
organized as follows. In Section 2 we study CTRW under a general waiting time pdf and 
two specific waiting time pdf's. In Section 3, we show the exact solutions for pdf p(x,t). 
Finally, conclusions are presented in Section 4. 

II. CONTINUOUS TIME RANDOM WALK AND INTEGRO-DIFFERENTIAL 
EQUATION FOR GENERAL WAITING TIME PROBABILITY DENSITY FUNC- 
TION 

Eq.([7j) was already employed to study diffusion behavior when the form of waiting time 
pdf g(t) was first specified such as exponential function or a power- law function in the long- 
time limit 4|. However, for a generic form of g(t) Eq.flTJ) is not convenient to be used to 
study diffusion behavior due to the difficulty in the inverse transformations of Fourier and 
Laplace in Eq.([7j). In this paper, we will derive an integro-differential equation for CTRW, 
which is valid for general waiting time probability density function. The details of derivation 
can be found in Appendix A. This integro-differential equation for CTRW is written by 

d P (x,t) [* n Wi), n d f* n 9MMi).. nm 

9{t-ti)—z: dt 1 = D— g(t-h) — dfr. (10) 



dt J dti dt Jq dx 2 

For the case of g(t) = t a ~ 1 /r(a), it is noted that the Caputo fractional derivative appears 
on the left side of Eq. ([TO]) and the Riemann-Liouville fractional derivative appears on the 
right side of Eq. ( TTOl) . The Caputo fractional derivative requires the integrability of the 
derivative and contains the initial value of the function. Therefore the Caputo fractional 
derivative is more restrictive than the Riemann-Liouville fractional derivative is. It should 
be pointed out that the use of these operators may lead to different behaviors including 



unphysical behavior in different systems 34|, |35[. Eq. (II 01) is derived from a well-defined 



physical process, and no ambiguities and unphysical behaviors are expected. The left side 
of Eq. ( II Op shows the variation of p(x, t) with respect to time, which does not depend only 
on the ordinary derivative operator but also the difference between ordinary and non-local 
integral operators. 

Our Eq. (10) is different from the fractional Fokker-Planck Eq. ([I]) in the structure 
because the fractional Fokker-Planck Eq. ([I]) does not have the second term of the left side 
of Eq. (HOD. Substituting g(t) = t a - l /T(a) into Eq. (jTHJ cannot directly obtain Eq. JIJ. In 
order to exclude the second term of the left side of Eq. (ITU1) . the asymptotic behavior for 
g(t) must be carefully treated. 

We note that 

X - dx . (11) 



dt J^ dt 

Substitute Eq. ((IDj) into (TIT]) yields 

d(x 2 ) f* , ,d(iV 3 f* , . f°° 2 d 2 p(x,t) , , 

After integrating the second term of the right side of Eq. (IT21) by parts twice, we have 

d(r 2 ) /"* d (r 2 ) 8 /"* Z" 00 

dt = / ^ (t " tl} dt dtl + 2D lk / S(t-*i) P& ^) d ^i . (13) 

where we also consider that linx^-i-oo p(x, t) and it decreases faster than 1/x. The normal- 
ization of the pdf p(x,t) requires that j_ p(x,t)dx = 1. Finally, Eq. (ll3j) can be rewritten 

as 



d{x 2 ) 



* jfsC-'O^ + ^ljf^t-Odt,. (14) 

Employing the Laplace transform in Eq. ( TH1) we obtain 

{X)L -— + Sil- 9 (S)} (15) 

in Laplace space. We have shown that from Eq. ( TTOl) even without the previous knowledge 
of the waiting time pdf, the second moment of displacement i.e. the diffusion behavior can 
also be obtained in Eq. (TT5|) . This is one of advantages of Eq. (TT0|) . Therefore, Eq. (TT5l) is 
very general and valid for any waiting time pdf. 

In fact, different diffusion behaviors can be found when substituting different waiting time 
pdf's into Eq.(TT5l). Considering that the waiting time pdf must be positive and normalized, 



there exist only a few simple functions that can be used as waiting time pdf's. Now we 
study diffusion behaviors in two specific cases of waiting time pdf's. 

First case: g\{t) = \e~ xt . 

In this case the characteristic waiting time is finite. The Laplace transform of gi(t) is 
given by 

9i(s) = t^- • (16) 



X + s 
Substituting Eq. (TT6|) into Eq. ([15]) we obtain 

(x 2 ) = (x 2 ) + 2D\t . 



(17) 



Eq. ffTTl) shows the normal diffusion behavior. It has also shown that the generalized diffusion 
equation ([TO]) can reduce to describe the normal diffusion process when the waiting time pdf 
is a simple exponential function. 

Second case: g2(t) = \ a t a ~ l E aya {— X a t a ). 



E^ u (y) is the generalized Mittag-Leffler function defined by 



36. 
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E», v {y) = J2 



If 



n=0 



r (z/ + fin) 



H > 0, v > . 



(18) 



We restrict the parameter a to < a < 1 because E a ^ a (—X a t a ), in this interval, is completely 
monotone for t > [36]. We note that #2 (t) can also be written as 

d 



92(t) 



dt 



E a A—X a t c 



This waiting time pdf has been used to link the CTRW model and fractional master equation 



39 



40J. Fulger et al. 41] incor 



381 ] . It has also been adopted as subordination function 
porated this waiting time pdf into Monte-carlo simulation of in uncoupled continuous-time 
random walks with a Levy a-stable distribution of jumps. The Levy a-stable distribution 



of jumps used in 



41] has a diverging jump length variance in contrast to the finite jump 



length variance given by Eq. ([6]) in this paper. 

The waiting time pdf g2{t) has the following asymptotic behavior 

1 



92(t) 



(19) 



\ a t 1+a T(-a) ' 

which has a power-law behavior. The characteristic waiting time is divergent. The Laplace 
transform of g^ (t) is given by 

92(s) = v^ . (20) 



Xa + S c 



Substituting Eq. (120]) into Eq. (fT5|) we obtain 

2D\ a 



(»>) = (» 2 ) + T^f- < 21 > 



I, 



This result is the same as that one obtained in [4j for long-time limit using approximation 
of a long-tailed waiting time pdf. However, it should be noted that Eq. (12"T|) describes the 
subdiffusive processes for the all time regions rather than only for long-time or short time 



i. 



region. This point is different from the case of long-tailed waiting time pdf in Ref. |4[]. 



III. EXACT SOLUTION FOR PROBABILITY DENSITY 

Generally speaking, using Eq.(j3]) we can calculate the pdf in the framework of CTRW 
model. However, it is hard to obtain the exact pdf from Eq.([3]). Now, we can obtain the 
exact solutions for p(x,t) for the two waiting time pdf's in Section 2 from Eq.([7]). In the 
following, we consider the initial condition as p${k) = 1. After we do Fourier inverse for 
p{k,s) in Eq.(J7]), we have 

1 _ a ( s ) r°° e ikx 
p(x,s)= J { /, / —^&k. (22) 

We note that the denominator can have three different poles: 1 — g(s) = gives trivial 
result, (1 — g(s))/g(s) < gives the poles on the real axis and (1 — g(s))/g(s) > gives the 
poles on the imaginary axis. For the waiting time pdf's given in the previous section we can 
restrict to the case of the poles on the imaginary axis. Thus, the solution for p(x, s) is given 

by 

1 l — q(s) __l£L/IEsgr , n 
p(x,s =_- _J_J^ e W im . 23 

2ttVDs V g{s) 
Before obtaining the exact solutions to p(x, t) we can calculate the second moment in Laplace 

space which yields 

2 2Dg(s) 

{X)l ~ s[l-g(s)]- [M) 

This result is similar to ffT5|) . except the term of initial value (a; 2 ) . This indirectly shows 
that Eq.(l23l is correct. 

For the case of gi{t) = \e~ xt the solution for p(x, t) is given by 

1 x 2 

yAnDXt 

7 



The pdf shows the Gaussian shape and it is exactly the solution of the normal diffusion 



equation 



\i. 



For the case of 02 (0 = ^at a l E a>a (—\ a t a ) the solution for p(x,t) is given by 

A _ 1 ^ (-i) n f x 2 y 



This result is exactly the solution of the fractional diffusion equation (P). The pdf 
presents pronounced cusps in x-coordinate for different times. This form is different from 



14 



the Gaussian shape (j25|) which presents smooth shape in the same coordinate |4|. 

Eq. ( 126]) can also be expressed in terms of the Wright function <& v j (y) J42| in the following 

form: 

1 ( \x\ \ 

p 2 (x,t) = $_^ !_« ' ' 5 , (27) 



where 



*«fe) = E^(|r^y. /5>-l.*eC. (28) 



n=0 



IV. CONCLUSION 

In this work we have investigated the CTRW model with the decoupled jump pdf. We 
have derived the integro-differential equation (TlO|) . From Eq. liTUj) we have derived diffusion 
behavior that is explicitly related to the waiting time pdf. Eq. ffTQj) is more flexible to be 
used than Eq.fl7J) because using Eq.(10) we can avoid Laplace and Fourier transformations 
and their inverse operations, for instance, Eq. (114p is ready to be used for calculating the 
second moment numerically. Both Eqs.f llOj) and (1151) are very general and valid for any 
waiting time pdf. To our knowledge, Eqs. (TT0|) and (TT5|) have not been published elsewhere. 
In fact, Eq. ([7|) or Eq. ffTUj) describes random walk model without external force. However, 
interesting physical problems involve external force. Generalization of Eq. (llOp to include 
the effect of external force is a natural step. In addition, possible generalization of Eq. ffTU]) . 
e.g., using the master equation approach [4j, to include the effect of external force would be 
more convenient and appropriate than Eq. (J7J). 

Using Eq. fflOl) we have studied diffusion behaviors for two cases: (1) in the case of 
gi(t) = \e~ xt , a normal diffusion or normal Brownian motion is found; (2) in the case 
of gi{t) = ^ a t a ~ l E at0t (— \ a t a ) we have discovered a subdiffusive behavior for the all time 



8 



regions rather than only for long-time or short time regions that appears in the literature. 
We have also obtained the exact solution for pdf's in the two cases. For the first case, 
the exact solution of pdf is Gaussian. For the second case, the exact solution of pdf is 
non-Gaussian, but it can reduce to Gaussian pdf when a, — 1. 
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Appendix A: Integro-differential equation 



In order to obtain the integro-differential equation (fTUj) from Eq. ([7]) we employ the 
Fourier and Laplace transforms of the following expressions: 

~d 2 p(x,t)' 



T 



dx 2 



-k 2 p(k,t) 



and the convolution theorem for Laplace transform 



C 



git-t^ p(x,ti)dt 1 



uo 



g(s)p(x,s). 



Firstly, we use the inverse Fourier transform and Eq. flAip in Eq. ([7]), and we find 

d 2 p(x, s) 



sp(x } s) - sg(s)p(x, s) - Dsg(s)- 



dx 2 



[l-g(s)]p(x,0). 



After both sides of Eq. (IA3P are divided by s, we reorganize it as follows: 

1 1 d p(x s) 

p(x, s) - -p(x, 0) - g{s)p(x, s) + -g(s)p(x, 0) = Dg(s) — —j— 
s s ox z 



(Al) 



(A2) 



(A3) 



(A4) 



Using Eq.( 1A2l) and applying the inverse Laplace transform in Eq.( 1A4l) we have 
p(x,t)-p(x,0)- [ ( ? (t-t 1 )p(x,t 1 )dt 1 +p(a;,0) / g{t 1 )fo x = D [ g(t-h) T^ dfr 



dx 2 



(A5) 



Using the Leibnitz theorem for integral, we have the following relationship: 



0_ 

dt 



g(t-ti)p(x,t 1 )dt 1 



g(0)p(x,t) + 



d_ 

Of 



ff(*-ti)p(x,ti)d*i. (A6) 



The integral on the right side of Eq. (1A6|) can be expressed as 



* d f* d 

g(t-ti)p{x,ti)dti = - / —-g(t-ti)p(x,ti)dt 1 . 



(A7) 



The integral on the right side of Eq. (lA7p can be integrated by parts, and Eq.( 1A7p can be 
further written as 

/—#(*- ti) p{x, h)dti = -g(0)p(x, t) + g(t)p(x, 0) + / g (t - h) —p{x, i 1 )dt 1 . (A8) 

Substituting Eq. (IA80 into Eq.( lA6p . Eq.( 1A6p can be rewritten as 



d_ 

dt 



g(t-ti)p(x,t 1 )dt 1 



g{t)p{x, 0) + J g {t - h) —p(x, t^dh. (A9) 



Firstly employing the operator d/dt to Eq.( lA5P and then substituting Eq.f ]A9P into the 
equation, finally we have 



dp(x,t) 
dt 



f* n 9^ti), n d f* n ^M, 



(A10) 
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